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Cloning and expression analysis of a UBX domain-containing 


protein gene LmUBX2 from the migratory locust, 


Locusta migratoria manilensis ( Orthoptera: Locustidae ) 
RAN Yong-Hong, HUANG Ting, HAO You-Jin, CHEN Bin, HE Zheng-Bo “ 


(Institute of Entomology and Molecular Biology, Chongqing Normal University, Chongqing 400047 , China) 


Abstract: [ Aim] Ubiquitin regulatory X (UBX) domain-containing proteins function as cofactors of 
p97/CDC48 which is involved in multiple ubiquitin-related processes including protein degradation 
through the ubiquitin-proteasome system and homotypic membrane fusion. Our study aims to clone the 
UBX domain-containing protein gene from Locusta migratoria manilensis, and to analyze its expression 
patterns in different tissues and at different developmental stages, so as to provide the foundation for 
further function research. [Methods] A transcriptome database for L. migratoria manilensis was mined 
through bioinformatic analysis. Gene expression levels were analyzed in different tissues of L. migratoria 
manilensis adults and at different developmental stages using real-time PCR. [Results] A UBX domain- 
containing protein gene was identified, and named LmUBX2. The open reading frame (ORF) encodes a 
protein of 339 amino acids with a predicted molecular weight of 37.8 kDa and a theoretical pI of 6.03. 
BLAST analysis showed that it shares 37% -64% amino acid sequence identities to other UBX domain- 
containing proteins. LmUBX2 possesses a conserved N-terminal UBA domain and a C-terminal UBX 
domain. Sequence comparison and phylogenetic analysis revealed that LmUBX2 is a member of the 
SAKS1 subfamily. Gene expression analysis results showed that LmUBX2 was expressed during the whole 
life cycle, and highly expressed in ovary and testis. [Conclusion] These findings suggest that LmUBX2 
might be involved in multiple physiological processes of L. migratoria manilensis. Especially, LmUBX2 
might be related to the reproduction of L. migratoria manilensis, although further research needs to be 
performed to confirm this correlation. 

Key words: Locusta migratoria manilensis; UBX domain-containing protein; LmUBX2; gene cloning; 
expression pattern 
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1 INTRODUCTION 


Ubiquitin regulatory X (UBX) domain was first 
described in the hypothetical human protein Y33K 
1996 ). Subsequent 


UBX domain-containing 


and Bucher, 
that 

proteins are widely distributed in all eukaryotes 
( Schuberth and Buchberger, 2008). Structural 
analysis of the UBX domain of FAF-1 and p47 reveal 
that they are structurally homologous to ubiquitin and 


( Hofmann 


researches showed 


act as a p97-binding module ( Buchberger et al., 
2001; Yuan et al., 2001; Dreveny et al., 2004). 
The p97, also known as VCP in higher eukaryotes 
and Cdc48p in yeast, is a homohexameric ATPase 
(ATPase associated with a variety of activities ) 
involved in multiple cellular processes, including 
protein degradation through the ubiquitin-proteasome 
system and homotypic membrane fusion ( Woodman, 
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2003; Wang et al., 2004; Ye, 2006; White and 
Lauring, 2007). 

UBX domain-containing proteins have been 
identified in a wide range of species (Schuberth and 
Buchberger, 2008 ). Based on the sequence 
homology outside the UBX domain, UBX domain- 
containing proteins can be divided into several 
conserved subfamilies, including UBXD-1, UBXD- 
3, FAF1, p47, SAKS1, TUG, Rep8 and so on 
(Schuberth and Buchberger, 2008). Caenorhabditis 
elegans possesses 
(UBXN-1 to -6) , which were dynamically expressed 


six UBX-containing proteins 


throughout worm’ s development and in different 
tissues tested ( Yamauchi et al., 2007). All six 
UBXN proteins are cofactors for CDC48/p97 and 


induce embryonic lethal and sterile by the 
simultaneous knockdown of wbx«n-1, ubxn-2 and 
ubxn-3 ( Sasagawa et al., 2010). Caspar, a 


homologue of human FAF1, was found to be a 


APERE H (31371992); ERTL AA SRT 5 “IT” HP Ls Be OE H 
KŻ, &, 1990 FE, HEPA, WEWRE, WERE TT Tay Be PPE Ae, E-mail; 1048765288 @ qq. com 


362 BB eR Acta Entomologica Sinica 58 4 





specific negative regulator in the immune deficiency 
(Imd) signaling pathway of Drosophila by inhibiting 
the Dredd-dependent cleavage of Relish in vivo (Kim 


et al., 2006 ). 
antimicrobial peptide ( AMP) gene diptericin was 


In the caspar mutant flies, the 


constitutively expressed even in the absence of 
infections. While, the overexpression of caspar was 
sufficient to suppress the Imd pathway but not the 
Toll pathway (Kim et al., 2006). The knockdown of 
Locusta migratoria FAF1 (ImFAF1) by RNAi 
AMP gene 
prolixicin , suggesting that it shares a similar function 
as Caspar in Drosophila (He et al., 2013a). In 


major 


caused the hyperactivation of the 


malaria vectors Anopheles stephensi and 
Anopheles albimanus, the RNAi-mediated silencing 
of caspar efficiently prevents the development of 
malaria parasite Plasmodium falciparum (Garver et 
al., 2009). 

The oriental migratory locust, Locusta 
migratoria manilensis ( Meyen), is an orthopteran 
pest and a 
hemimetabolous insects for biological studies ( Kang 
et al., 2004; Ma et al., 2006). In previous studies, 


we identified 


representative member of the 


several UBX  domain-containing 
proteins from L. migratoria manilensis including 
LmFAF1 and LmUBX1, which were classified into 
the FAF1 subfamily and the UBXD1 subfamily, 
respectively, based on the sequence characteristics 


( He et al., 2013a, 


However, there is no more information 


and phylogenetic analysis 
2013b ). 
about other subfamilies in L. migratoria manilensis. 
In this study, a gene encoding the UBX domain- 
containing protein of the SAKS1 subfamily was 
identified through mining the transcriptome database 
(unpublished data) and designated as LmUBX2. Its 
expression patterns were further investigated at 
different developmental stages or in different tissues 
using real-time RT-PCR. Taken together, these 
results provided an insight into the possible roles of 
LmUBX2 in locust development. 


2 MATERIALS AND METHODS 


2.1 Experimental insect 

Locusta migratoria manilensis ( Orthoptera: 
Acrididae ) was originally collected from Hebei 
Province, China, and reared as previously described 
(He et al., 2006). 
2.2 RNA extraction 

Total RNA was extracted from different tissues 
( TaKaRa Dalian, China ) 
according to the manufacturer’ s protocol and treated 
with On-Column DNase to DNA 


contamination. RNA concentration and quality were 


using Trizol reagent 


remove 


estimated by spectrophotometry and agarose gel 
electrophoresis , respectively. cDNA was synthesized 
from 1 ug mRNA using the Smart Race cDNA 
Amplification kit (Clontech, Mountain View, CA, 
USA) according to the manufacturer’ s instructions. 
2.3 Cloning and bioinformatic analysis of 
LmUBX2 

The transcriptome database of L. migratoria 
manilensis was analyzed using human SAKSI 
(Y¥33K; LOC51035 ) as a query. The retrieved 
sequences were further analyzed by BLASTX against 
the non-redundant database at the NCBI. The open 
reading frame ( ORF) was predicted by ORF Finder 
(http://www. ncbi. nlm. nih. gov/projects/gorf/ ) . 
The theoretical isoelectric point and molecular weight 
were predicted using ProtParam and Compute pI/Mw 
(http ://web. expasy. org/compute_pi/). The signal 
peptides were predicted by the SignalP 4. 0 server 
( http://www. cbs. dtu. dk/services/SignalP/ ). 
Motif identification was carried out by ScanProsite 
(http ://prosite. expasy. org/) (De Castro et al., 
2006). Sequence similarity was analyzed by BLAST 
searching in GenBank (http://www. ncbi. nlm. nih. 
gov/). Multi-sequence alignment was generated 
using Clustal X (Thompson et al., 1997). Homology 
modeling was performed by the SWISS-MODEL 
server (http://swissmodel. expasy. org/) ( Biasini 
et al., 2014). An unrooted phylogenetic tree was 
constructed by MEGA 4.0 (Tamura et al., 2007 ) 
using the neighbor-joining method with 1 000 
bootstrap replications. 
2.4 Expression patterns of LmUBX2 at different 
developmental stages and in different tissues 

Gene expression was analyzed at different 
developmental stages including embryo ( mixed with 
various developmental stages), Ist, 2nd, 3rd, 4th 
and 5th instar nymphs, and adult. Locusts were 
collected in whole bodies (remains in the alimentary 
canal were removed) at 2 d after moulting, and 5 to 
10 individuals were collected for each sample. 
Meanwhile, gene expression was also examined in 
the midgut, ovary, testis, integument and fat body. 
Ovary and testis were dissected from 2 day-old female 
and male adults. Other tissues were collected from 
female adults. The locusts were dissected on ice, and 
the tissues were isolated and stored at — 80°C. 

Real-time PCR was carried out using SYBR 
PrimeScript RT-PCR Kit (TaKaRa, Dalian, China) 
and CFX Connect™ (Bio-Rad). First-strand cDNA 
was synthesized from 1 wg of total RNA using 
random hexamers. Gene specific primers were listed 
in Table 1. The PCR conditions were as follows; 


95°C for 15 s; 40 cycles of 95°C for 10 s, 58°C 15 s, 
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72° 20 s; and then 95°C for 15 s and 60°C for 30 
s. The Armadillo was simultaneously quantified as an 


internal control. The relative expression level was 


measured using the 27“ method ( Livak and 
Schmittgen , 2001). 


and the experiment was repeated three times. 


Each sample had three replicates 


Table 1 Primers used in this study 


Primer name 


Sequence (5' —3’) 


Purpose 





UBX2-F CTGCCTCAACTATGCCTACAC 

UBX2-R GCCTAACTGCTGCCAACTG 
18SF TTAAGCCATGCATGTCTAAGTA 
18SR 


LmUBX2 real-time PCR 
LmUBX2 real-time PCR 
18S rRNA real-time PCR 
185 rRNA real-time PCR 


TCTCAGGCTCCCTCTCCGGAATCG val-ti 


3 RESULTS 


3.1 Cloning and characterization of LmUBX2 
A gene fragment ( Unigene38793 ) was retrieved 
and showed 40% - 64% identity to other UBX 
domain-containing proteins at the amino acid level. 
This gene was named as LmUBX2 and was deposited 
into GenBank (accession no. KP280173). LmUBX2 
contains a 1 020-bp ORF encoding a 339-aa protein 
with a theoretical MW of 37.80 kDa and a pl of 6. 03. 
The genomic sequences of LmUBX2 were identified by 
BLAST against the L. migratoria whole genome 
shotgun (WGS) sequences ( Wang et al., 2014). The 
entire ORF of LmUBXZ2 is distributed in 6 contigs and 
consists of 7 exons (Fig. 1; A). Except the intron 1, 
the length of other 5 introns remains unknown due to 
the discontinuous genomic sequences (Fig. 1; A). 
Secondary structural prediction showed that 


A contig_ contig —-contig_ 

a f 360103 360105 360106 
10 001 b GAD 8239 

LmUBX? --- ME -----: N- -: -S -+ -+ 


B (56%-79% identity (51%-— 
to ae iia 


83% identity C a 
to other UBX) 


LmUBX2 possess 51. 33% 
extended strand, 3. 54% 
Therefore, 
represent the dominant secondary structures. Domain 
that LmUBX2 has 
functional domains including an N-terminal 
ubiquitin-associated ( UBA ) domain (aa 2 - 40), 
and a C-terminal UBX domain (aa 253 - 336), 
which respectively exhibited 56% - 79% and 
51% -83% sequence identities to other UBA and 
UBX domains (Fig. 1; B). The modeling structures 
of UBA and UBX domains of LmUBX2 and human 
SAKS1 were constructed through on-line SWISS- 
MODEL service (Fig. 1: E, F, G and H). 
Comparative analysis showed that the structures of 
UBA and UBX domains are highly conserved ( Fig. 
1), which are also consistent with the conserved 
amino acid sequences of UBA and UBX domain 


(Fig. 1: C, D). 


a-helix, 5. 90% 
B-turn and 39. 23% 
random coil. a-helix and random coil 


analysis revealed multiple 


contig _ contig_ contig 

360111 360112 360114 

me 783 a (14 895 i (10 115 b a 
Pa 241 S iss 7 48 


T. ODVN' 
ELTALES i a 


UBX 
mumo sD “PFE PAER TEER: 
GEE 
BA UBX 


Hssaksı MET 31200 


E TNES LVPSAVVI 
ak $ MERP = MAARLETR 


Ne yee E 


Fig. 1 


Genomic structure and domains analysis of LmUBX2 of Locusta migratoria manilensis 


A; Genomic structure of LmUBX2. The bold dotted lines represent the gaps between the contigs, while the thin dotted lines represent introns, or the 5- 


terminal and 3-terminal untranslated regions. “ 


the length of the contigs, while the figures below the bold grey line represent the length of the exons. 


domain of LmUBX2 and human SAKS1. 


D: Alignment of UBX domain of LmUBX2 and human SAKS1. 


?” indicates that the length of the intron remains unknown. The figures above the bold black line represent 


B: Domains of LmUBX2. C; Alignment of UBA 


Identical residues ( > 50% ) are shaded with 


dark gray. The comparison of the modeling structures of UBA domain LmUBX2 (E) and SAKS1 (F). The comparison of the modeling structures of UBX 


domain of LmUBX2 (G) and SAKS1 (H). LmUBX2; Locusta migratoria manilensis , 


Ubiquitin regulatory X domain; UBA ; Ubiquitin-associated domain. 


KP280173; HsSAKS1; LOC51035; UBX; 


Homo sapiens, 
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BLAST analysis showed that LmUBX2 is a 
relatively conserved protein sharing 37% -64% amino 
acid sequence identities with other UBX domain- 
containing proteins from plants, worms, insects, 
mammals and fungi. Amongst them, LmUBX2 is most 
closely related to the homologue (GenBank accession 
no. XP_005814923) from Xiphophorus maculatus with 
64% amino acid sequence identity. Phylogenetic 
analysis of the 17 UBX domain-containing proteins, 
based on amino acid sequences, is presented in Fig. 
2. The tree topology indicated that LmUBX2 was 
clustered with other SAKS1 homologues, suggesting 
that LmUBX2 is a member of the SAKS1 subfamily. 
While LmUBX1 andLmFAF1 were grouped into the 
UBXD1 FAF1 


respectively (Fig. 2). 


subfamily and the subfamily , 


100; LmUBX2 
92 ~ ZnUBX1 
98| |_______ CeUBXN-1 
53 |__ pss aKs1 
1A SpYDFB 
HsErasin 
54 CeUBXN-4 
CeUBXN-6 
100) m LmUBXI 
24| 93 HsUBXD1 
es FY. 
82 c HsFAF1 
HsETEA 
HsRep-8 | Rep-8 subfamily 
Hsp47 | p47 subfamily 
64 HsTUG | TUG subfamily 
HsUBXD3 | UBXD3 subfamily 
Fig. 2 Phylogenetic tree of LmUBX2 of Locusta migratoria 




















SAKS1 subfamily 





























UBXD1 subfamily 

















FAF1 subfamily 























76 


























manilensis and homologues of other species constructed 

by MEGA 4 with the neighbor-joining method based 
on amino acid sequences 
The number at each node indicates the percentage of bootstrapping after 
1 000 replications. LmUBX2; Locusta migratoria manilensis UBX2, 
KP280173; ZnUBX1: Zootermopsis nevadensis UBX1, KDR19499; 
CeUBXN1; Caenorhabditis elegans UBXN-1, CCD69888; HsSAKS1: 
Homo sapiens SAKS1 , LOC51035; SpYDFB; Schizosaccharomyces pombe 
YDFB, SPACI7C9. 11c; erasin; H. sapiens, ABM90426; CeUBXN4; 
C. elegans UBXN4, NP_498856; CeUBXN6:; C. elegans UBXN-6, NP 
_500648. 2; LmUBX1; L. migratoria manilensis UBX1, JN661173; 
HsUBXD1: H. sapiens UBXD1, AF272893; LmFAFI:; Locusta 
migratoria manilensis FAF1, AB758246; HsFAF1: H. sapiens FAF1, 
QO9UNNS; HsETEA; H. sapiens, AAH14001; HsRep-8; H. sapiens 
Rep-8 , D83767; Hsp47: H. sapiens p47, BC002801; HsTUG; H. 
sapiens ASPL, AF324219. 1; HsUBXD3: H. sapiens UBXD3, P_ 
689589. 1. 


3. 2 Expression of LmUBX2 at different 
developmental stages 

LmUBX2 expression pattern was analyzed 
throughout the life cycle by real-time PCR. The 
highest expression level was found at adult stage 
(F =48. 1560, P =0. 0001), while the transcript 
levels did not significantly vary during the embryo 
and nymphal development ( F = 2. 0630, P = 


0.1412) (Fig. 3; A). 
3.3 Expression of LmUBX2 in different tissues 
The transcript levels of LmUBX2 significantly 
varied in different tissues ( F = 554. 435, P = 
0.0001) (Fig. 3: B). It was highly expressed in 
ovary, followed by testes. Besides, the difference 
between ovary and testes (t = 1.6949, P =0. 1653 ) 
was significant (t = 5. 8330, P = 0. 0043). The 
transcript level in the ovaries and testes were 6.77-, 
and 5. 37-fold higher than that in the integument, 
respectively. 
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IG MG FB TE OV 
Tissues 
Fig. 3 Gene expression profiles of LmUBX2 throughout the 
life cycle (A) and in different tissues (B) of 
Locusta migratoria manilensis 
Armadillo was used as an endogenous control. EM; Embryo; N1; 1st 


instar nymph; N2; 2nd instar nymph; N3: 3rd instar nymph; N4; 4th 
instar nymph; N5; 5th instar nymph; AD: Adult; MG: Midgut; OV: 
Ovary; TE; Testis; IG; Integument; FB; Fatty body. Columns topped 
by different letters indicate significantly different means at the 0.01 level 
by ANOVA analysis. 


4 DISCUSSION 


In this study, LmUBX2 was identified and 
characterized from L. migratoria 
transcriptome dataset. Phylogenetic analysis 
indicated that LmUBX2 was clustered with other 
SAKS1 proteins, suggesting that LmUBX2 could be a 
member of the SAKS1 subfamily (Fig. 2). Previous 
studies thaliana , 
Saccharomyces cerevisiae, Homo sapiens, C. elegans 
genome contains 15, 7, 13 and 10 UBX domain- 
containing protein genes, respectively ( Rancour et 
al., 2004; Schuberth et al., 2004; Yamauchi et al., 
2007; Schuberth 2008 ). 


manilensis 


Arabidopsis 


revealed that 


and Buchberger, 
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Therefore, the three UBX  domain-containing 
LmUBX1, LmUBX2 and LmFAFI, 
obtained from the transcriptome data might not 


proteins , 


represent the entire repertoire of the UBX domain- 
containing protein family of ZL. migratoria 
Additionally, the UBX domain- 
containing protein genes of L. migratoria manilensis 
Fortunately, the L. 
migratoria manilensis genome has been sequenced 
and released ( Wang et al., 2014), which will 
facilitate the genome-wide characterization of all 


manilensis. 


remain to be discovered. 


UBX domain-containing proteins. 
Expression analysis results showed that 
LmUBX2 was throughout the 


developmental process of L. migratoria manilensis 


expressed 


(Fig. 3), suggesting that it might be involved in 
multiple physiological processes of L. migratoria 
LmUBX1 and ImFAF1 were also 
distinctly expressed at each developmental stage ( He 
et al., 2013a, 2013b). Of particular concern is that 
LmUBX2 was highly expressed in ovary and testis 
expressed in ovary and testis ( He et al., 2013b), 


manilensis. 


Similarly, LmUBX1 was mainly 


and LmFAF1 was simultaneously highly expressed in 
fat body, testis and ovary (He et al., 2013a). 
Taken together, our results suggested that these three 
UBX genes LmUBX1 , LmUBX2 and LmFAF1, are 
regulated in strikingly different manners. 

The function of UBX proteins of the SAKAI 
subfamily remains unknown. The C. elegans UBXN- 
1, the same subfamily with LmUBX2 , was expressed 
ubiquitously through all developmental stages 
(Yamauchi et al., 2007). Knockdown of ubxn-1, 
ubxn-2 and ubxn-3 induced embryonic lethal and 
sterile phenotypes ( Sasagawa et al., 2010 ). 
Furthermore, TRA-1A, a key factor in the sex 
determination pathway, accumulated in worms and 
inhibited spermatogenesis when ubxn-1, ubxn-2 and 
ubxn-3 were simultaneously knocked down. This 
suggests that UBXN-1, UBXN-2 and UBXN-3 may 
control spermatogenesis via the degradation of TRA- 
1A (Sasagawa et al., 2010). The higher abundance 
of LmUBX2 in ovary and testis indicates that 
LmUBX2 might be involved in the reproduction of L. 
further 
required to confirm these correlations. 


migratoria manilensis. But research is 


UBX domain-containing proteins function as 
cofactors of p97/CDC48, which is involved in 
multiple ubiquitin-related processes including protein 
degradation through the 
system. LmUBX2, like C. 
human SAKAI, 
domains, including an N-terminal UBA and a C- 


ubiquitin-proteasome 
elegans UBXN-1 and 


contains multi-ubiquitin-related 


terminal UBX domain (Fig. 1; B). UBX domain is 
structurally homologous to ubiquitin and acts as a 
whereas the UBA domain 
seems to be responsible for recruiting and stabilizing 
ubiquitinated proteasomal substrates (Schuberth and 
Buchberger, 2008). Existing results reveal that 
UBXN-1, UBXN-2 and UBXN-3 are likely involved 
in the degradation of TRA-1A (Sasagawa et al., 
2010). Moreover, ovary and testis cells are renewed 
relatively faster rate 
exuberant protein metabolization. Therefore, we 
speculated that LmUBX2 may be involved in the 
degradation of some proteins in ovary and testis of L. 
migratoria manilensis , but further studies need to be 
performed to test this hypothesis. 


p97-binding module, 


at a accompanying with 
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AE Kee UBX AMMEASSEARA 
LmUBX2 BY 5c BE FH Fe IK FP AT 
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(HUTA BE EE FET, EK 400047 ) 

















#32: [ AMY) UBX Ai yh FA p97/CDC48 Hyhi 














EH, WRZ R-E A AS R E I EJ eT] a SS ES Ee MER KA Locusta migratoria 








WAF. p97 AIZ REIRA & PP A Md E Pike a BS) 




















manilensis ( Meyen) ft) UBX AREA E AEN , PITRE AAC RAHA , Py EAE EE UBX 2 Pa a 9 EA 
FEA Be ee AE. [IA La TAR CE oR 

































































HCH oe Me UBX RE A E ETAL, OR AS INT et 






































PCR ERITA KES Te] AC A AS ER PIRA. [ER ] SEES KEA — A UBX 2 fey Sak 
LSE ASE , MAA LmUBX2. LmUBX2 FF BC BEHELK 1 020 bp, 4th 399 PAIR , EOL OP T EEEE DS 



































37.8 kDa #1 6.03 , 5th UBX REA E AY AER BE 37% ~ 64% ,N wig C SHI 























AESF HY UBA 


























ZAIRA UBX Sith. FI RAR RRA PT E LmUBX2 ja F SAKSI WAG. EETA ,LmUBX2 FEE 




















MEMIPA RIS (ELS HO Bea OP ic ; ETA R PRA 
Ho [WWR RH LmUBX2 HY EB AR KELPE ME, JHA AR 9 A 



































KIS , (ELLE AG ER AN ER PRAF te 
K AAA m IR 
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KHR: AKA KA; UBX MARUS EAH; LmUBX2; AE; RA 
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